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1 INTRODUCTION 



The origin of the intraday variabihty (IDV) (Wagner & 
Witzel 1995) seen in many flat-spectrum radio sources is a 
major astrophysical puzzle. Such rapid variability, if intrin- 
sic to the source, implies a very small source linear extent 
and correspondingly extreme values for the source surface 
brightness temper ature, Tb. For exam ple, in the BL Lac 
object 0716+714 (Wagner et al. 1996), where an intrinsic 



origin is favoured due to the correlation between the opti- 
cal and radio variations, the implied brightness temperature 



ABSTRACT 

Radio monitoring of the gravitational lens system B0218+357 reveals it to be a 
highly variable source with variations on timescales of a few days correlated in both 
images. This shows that the variability is intrinsic to the background lensed source and 
suggests that similar variations in other intraday variable sources can also be intrinsic 
in origin. 
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small (of order ^arcsec), the source brightness temperature 
can still greatly exceed the Compton limit, by several orders 
of magnitude in the case of PKS 0405-385. 

Even though both intrinsic and extrinsic explana- 
tions have their difficulties, it would represent considerable 
progress if it could be proven that the variations had one or 
the other origin. The study of radio sources which have mul- 
tiple gravitational images offers an elegant way to do this. If 
an IDV source is multiply imaged, only intrinsic variations 
should be correlated in the images. 



is of o iidei 10 K, well in exLess of the 

limit or • ■ lO^'^K for incoherent oynohrotron oourcoa. Invok 



The gravitational lens system B0218+357 (Patnaik et 



al. 1993) has been subject to many monitoring campaigns 



ing relativistic boosting can reduce the implied brightness 
temperatures, but often the Doppler factors req uired are so 



high (^100) as to be theoretically uncomfortable (Begelman 
Rees & Sikora 1994) and inconsistent with the s tatistics of 
superluminal motion (Vermeulen & Cohen 1994). 



Extrinsic effects that can cause rapid flu x density varl 
ations include interstellar scinti llation (ISS) (iRickett 199( ) 
and gravit ational microlensing ( ISopal-Krishna^T'Subrania- 



with the ultimate goal of measuring the Hubble parameter 
using the method of Refsdal (1964). The system consists 
of two images (A and B separated by 335 mas) of a com- 
pact flat-spectrum radio source and a steep-spectrum Ein- 
stein ring (Fig . [l|). The redsh ift of the background lensed 
object is 0.96 (Lawrence 1996). The flux densities of A and 



nian 1991). Whilst the latter has not been suggested as 



a major contributor to the IDV of any source to date, 
ISS has recently been successful in explaining the rapid 
variability of several sources, perhaps most notably th e 
quasars J1819-I-384 5 ([Dennett-Thorpe fc de Bruyn 2"00c| ) 
and PKS 0405-385 ( [Kedziora-Chudczer et al. 199?! ). Where 
ISS is the dominant cause of the variations the source bright- 
ness temperature can be reduced, to that requiring a reason- 
able Doppler factor (~10) in the case of J1819+3845. How- 
ever, as the source size required for scintillation is itself very 



B are approximately 800 mjy and 200 mjy respectively and 
typically have ~10 per cent linear polarizations. We have 
recently measured a time delay of 10.5 z fc 0.4 d in this sys - 
tem exploiting its rapid radio variability (Biggs et al. 1999). 



In this note we use our existing monitoring data and show 
that variability with an implied brightness temperature of 
^lO^^K is correlated in the two images. We argue that this 
cannot be due to ISS in the Galaxy or microlensing in the 
lensing galaxy and must therefore be intrinsic to the back- 
ground source. 
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0218+357 VLA 14.915 GHz 18/11/92 
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Figure 1. VLA 15-GHz radio map of B0218+357. As well as the 
two compact components (A to the right) and the Einstein ring, 
also clearly visible is a (non-lensed) radio jet to the south. 



2 OBSERVATIONS AND LIGHT CURVES 

B0218+357 was observed with the VLA in 'A' configura- 
tion between the months of 1996 October and 1997 Jan- 
uary. Observations were taken at two frequencies, 15 GHz 
and 8.4 GHz. In all, data were obtained at 47 epochs, with 
an average spacing between observations of ~2 d. Full de- 
tails of the observations and the data reduction are given in 
Biggs et al. (1999). Although the light curves are similar at 
both frequencies, it is at 15 GHz that the shortest timescale 
variability is most pronounced. Therefore we concentrate on 
the higher-frequency data in this paper and in Fig. ^ we 
reproduce the 15-GHz light curves for the two images. The 
time delay of 10.5 ± 0.4 d was determined from these and 
other light curves using a chi-squared minimisation tech- 
nique along with a Monte Carlo analysis of simulated light 
curves to derive the associated error. Knowing the delay 
and the relative magnifications of the images, it is possible 
to remove the time delay from the light curve of one of the 
images and produce a combined light curve. This is shown 
in Fig. ^. It is clear from a visual inspection of the combined 
light curve that the variations in the two images correlate 
very well. More formally, the time delay analysis of Biggs et 
al. (1999) showed that the reduced chi-squared (x^) of the 
fit was 0.7, indicating an excellent fit. 

We have no direct evidence for variability from the mon- 
itoring data on time-scales less than the typical sampling 
interval of ~2 days. However, in addition to the VLA moni- 
toring data we have a VLA (also A' configuration) 15-GHz 
10 hr dataset of B0218-I-357 from which the map in Fig. | 
was made. Observations of B0218-I-357 were taken in scans 
of approximately 20 minutes duration, each separated by a 
~2 minute scan of the phase-calibration source B0234+285. 
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Figure 3. VLA 15-GHz total flux density light curve constructed 
from combining the component A (open squares) and component 
B (filled squares) data shifted by a time delay of 10.5 d and scaled 
by the flux density ratio (A/B) of 3.73. 



This source was also used to calibrate the instrumental po- 
larization residuals. 

In order to study the short-term variability properties 
of B0218-F357 we have separated the data into the individ- 
ual scans (of which there are 20) before splitting each of 
these further into four equal sub-scans. This gives a total of 
80 epochs of da ta which were then model-fitted in difmap 
( Shepherd 1997 ) and in MPS using the uvfit task. The lat- 
ter was used for fitting to the Stokes Q and U parameters as 
well as to check on the difmap Stokes I results. The model- 
fitting was simplified by filtering the data in the (u, v) plane 
so as to remove the short baselines that respond to the dif- 
fuse Einstein ring emission. This Fourier filtering technique 
enabled us to restrict the model to two point-like compo- 
nents and a good fit (x'^ ~ 1) obtained. 

In order to avoid as much as possible errors in the flux 
calibration we have only considered the flux density ratio 
between A and B, the percentage polarization of each im- 
age and the difference in the polarization position angle be- 
tween the two components. These quantities should all be 
relatively free from systematic calibration errors. The un- 
certainty in our results has been estimated by calculating a 
mean and a standard deviation from the four epochs that 
make up each scan. We do not show any of the 'light curves' 
derived from this process as we do not believe that they 
show any evidence for significant variability (< 1 per cent) 
over the 10 hr period. 



3 DISCUSSION 

Does B0218+357 qualify as an intraday variable radio 
source? The variability in this lens syst em is less dramatic 
than that seen in the sources 0716+714 (Wagner et al. 199C) 



and 0917-1-624 (Rickett et al. 199E). An even greater contrast 



exists between B0218-(-357 and the quasars PKS 0405-385 
(tKedziora-Chudczer et al. 1997| ) and J1819+3845 ( [Dennett 
Thorpe fc de Bruyn 200(]| ) where much greater flux den- 
sity changes (over 300 per cent in J1819-I-3845) occur on 
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Figure 2. VLA 15-GHz total flux density light curves. Component A - left, component B - riglit. 



timescales of hours. However, the flux density of B0218+357 
does vary significantly over the typical time-scale of our 
monitoring observations, ^2 d, which is equivalent to ~1 d 
in the rest-frame of the source. In particular we draw atten- 
tion to the rapid fall in the A light curve that occurs towards 
the end of the monitoring where a decrease of about 60 mjy 
in the total flux density occurs over a period of around 5 d. 
This is associated with a rapid change in percentage polar- 
ization of about 1 per cent. Unfortunately, both for the time 
delay measurement and the present discussion, the event is 
too close to the end of the monitoring period for us to be 
able to see the corresponding feature in the B curve. Events 
that are obviously correlated in the light curves of Fig. ^ 
occur i n between days 20 and 40, the biggest of which is a 



transforming to a frame comoving with the source. Distance 
measurements were calculated assuming a Hubble parame- 
ter of 65 kms^^ Mpc~^ and a flat universe with Q, = 1. A cor- 
rection has also been made for the lens magnification using 
the model of Biggs et al. (1999). Although not stated in that 
paper explicitly, the absolute magnifications of components 
A and B are equal to 2.1 and 0.6 respectively (L.V.E. Koop- 
mans, private communication). We find values of Tb in the 
range lO^* — lO^^K. Brightness temperatures as high as these 
are greatly in excess of the inverse Compton limit of approxi- 
mately lO^'^K for incoherent synchrotron radiation. This can 
be reconciled by assuming bulk relativistic motion within a 
radio jet of Doppler factor D. Variability-derived brightness 



temp eratures are then reduced by a factor D (Readhead 



change of about 40 mJy over 10 d. 



1994) which for B0218+357 would require a Doppler factor 



From the observed variations it is possible to estimate 



of <10. This is consistent with observations of superlumi- 



minimn m brightncgg temperatures for the varying part of 
the ba ckground sourco. For the purposoE of thic diccuB — 



nal motion observed in the general source population (Vor 



meulen & Cohen 1994), but regular VLBI observations have 



sion we shall restrict ourselves to consideration of the most 
prominent and defined feature in the light curves, the change 
in fiux density, AS, of 60 mJy near the end of the moni- 
toring. Due to the uncertainty in defining the timescale of 
flux density variations, t, we have estimated this parame- 
ter in several different ways. Firstly, and most simply, we 
have taken t to be equal to the duration of the feature. 
At, as measured directly from the light curves i.e. approx- 
imately 5 d. Another approach is to weight At by the flux 
density of the varying part of the background source, 5vary 
i.e. t = At ISvary/AS*! (e.g. Wagner & Witzel 1995). What 
value should we assign to 5'vary? The conservative approach 
in terms of brightness temperature estimation would be to 
set this equal to the flux density of the component, ~ 1 Jy 
i n the case of component A. Ho wever. VLBI observations 
(Patnaik, Porcas & Browne 1995) have shown that approxi- 
mately one third of the flux density of components A and B 
originates from a jet-like subcomponent that is unlikely to 
vary signiflcantly over the course of our monitoring observa- 
tions. Therefore, Svary and t are reduced accordingly. 

We have converted our variability timescales into 
brightness temperatures using standard formulae such as 
equation (8) of Gopal-Krishna & Subramanian (1991), 



yet to reliably detect any change in the core-jet substructure 
of the A and B images of B0218-I-357 (e.g. Porcas & Patnaik 
1996). 

Could the correlated variability seen in the images of 
B0218+357 be produced by any of the extrinsic mecha- 
nisms? First let us consider interstellar scintillation. Let us 
assume that a typical velocity of the local Galactic inter- 
stellar medium with respect to the Earth is 20 km s~^. To 
produce correlated changes in images A and B with a time 
lag of <10 d the screen would have to move >335 mas in 
this time. With the above velocity the scintillating screen 
would have to be at a distance of less than 0.3 pc in order 
to produce correlated variations. This then allows us to dis- 
count the ISS explanation in two ways, the flrst of which is 
to note that 0.3 pc is very much less than the scale height 
of the Galactic disk. Secondly, at this distance and at this 
frequency the angular size of the flrst Fresnel zone is 0.1 mas 
(e.g. Walker 1998). Any source smaller than this would nec- 
essarily scintillate, but this has not been observed (VSOP 
calibration sources for example; J. E. J. Lovell, private com- 
munication). Walker (1996) has also concluded that as the 
typical separations (^1 arcsec) of galactic-scale mass gravi- 
tational lens systems are large compared to the typical size 
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of interstellar irregularities, any scintillation of the multiple 
images will be independent. Although the microlensing ex- 
planation might appear tempting, due to the very obvious 
presence of a lensing galaxy between us and the background 
source, the lines of sight when passing through this galaxy 
are so far separated that a correlated time delay of ~10 d can 
be completely ruled out. Microlensing in the host galaxy of 
the background radio source can be discounted for the same 
reason as scintillation in the Galaxy. 



4 CONCLUSIONS 

The source in the lens system B0218-I-357 is bright and 
rapidly variable with variations on a timescale of a few days 
that are correlated in both images. Thus on these timescales, 
at least, the mechanism producing the major variations must 
be intrinsic and any contribution from interstellar scintilla- 
tion small. Therefore, although recent results have favoured 
an extrinsic origin for the intraday variations seen in some 
sources (particularly J1819-)-3845 and PKS 0405-385), our 
results would suggest that IDV in these and other objects 
can also be partly intri nsic. For exam ple, 230 GHz observa- 
tions of PKS 0405-385 (Wagner 1995) have shown rapid vari- 
ability (with an associated brightness temperature of lO^'^K) 



that cannot be due to scintillation. Clearly more intensive 
monitoring of B0218-I-357, possibly in bursts separated by 
the time delay, would yield even tighter constraints on the 
intraday variability and enable the time delay to be refined. 
The same argument used here could also be applied to other 
variable lens systems, for exa mple PKS 1830-211 (Jauncey 
et al. 199l|; |LoveU et al. 1998|). 
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